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ABSTRACT 


Several techniques for finding approximate solutions to 
certain classes of nonlinear differential equations are in- 
vestigated. The nonlinear systems evaluated are second order 
HPeeaeauddrarte and Cubic terms, having driving functions and 
Pibelal COnditions. Pipes’ technique iS used to reduce the 
nonlinear differential equation to a system of linear dif- 
Pa Metalecoudtions. The Brady-Baycura technigue makes use 
Of nonlinear Laplace transforms to Soeain the solution. hve 
solutions are compared to the well known Runge-Kutta numerical 
ieewod SOlLUELON USing the digital computer. Greater accuracy 
was found uSing the Brady-Baycura method, but the simplicity 


eco ie miodmMakes Lt More attractive to the engineer. 





el 


ab. 


WABEE OF CONTENTS 


ie DOOM. 5 6 tlw lw lll ltl tll lll lt 
METHODS CIT ss es CC Ct 
A. PIPES' METHOD CUS ECCU—e 2. “ces. 6) 4) se 
Beemer —bpeaveURA METHOD .... « © © « « »« 
See lo U-OnTsCH-PIPES-NOWACKI METHOD .... 
PUSEMIRG@EUPEONG NS) 5 . ss kw ke wt tl 
PIC TIENEN (SPRING . . « «© «© «© © «© © © «© «@ 
1. Solution by Pipes Reversion Method . 
2. Solution by Brady-Baycura Transform 
RISTRWIGICMMEEEECEE TS .: ~@ 5 6 4 6 © «6 « © «© * 
3. Solution by Tou-Doetsch-Pipes-Nowacki 
DEIOUNGGIMECNI TNS) << 6 5 « sos © @ « * « 
4. Comparing the Methods oe ey 
B. NONLINEAR SPRING WITH DAMPENING Seed Somers 
iWeeoeluraon by Pipes Method .... . « 
ZoOolielom Oy Brady—-Baycura Method . . 
Saeecomparang the Solutions et 
C. NONLINEAR SPRING WITH CUBIC TERM ... . 
imoCOohMitron by Pipes Technique eens 
ZC Olutton by Brady—Baycura Method .. 
SweCoMpaciscon ©f Solutions See c cee cc wee 
D. NONLINEAR SPRING WITH DRIVING FUNCTION . 
mc Oltcrtenby Pipes Method .7. . .« . « 
Ze oem renhwoy Brady-Baycura Technique 
Scouparing the Solutions a ee 
E. NONLINEAR SPRING WITH CUBIC TERM AND 


FORCING FUNCTION eas sass 6 ek eC 


10 
10 


10 


13 


5 
ey 
by 
24 
Ze 
30 
30 
38 
42 
43 
4h 
Sal. 
a2 


56 


56 





os 


Ze 


St 


Solution by Pipes Method 


Solution by Brady-Baycura Technique 


Comparing the Solutions oars 


Ny CONCLUSIONS > CS. - SS = ree 


SeeeNDix A = 


See NDI xX B= 


meee NDiX C — 


BIBLIOGRAPHY 


hivemse Laplace Transform of 


Evaluation of A, in Forced 
Quadratic System >) Se 


Solution of Last Term of F. va 
Brady-Baycura Method = a : 


mye rtAGL DESTRIBUTION LIST SeMe no) 5: -s me Je 


FORM DD 


1473 


62 
63 
66 
72 


73 


74 


76 
78 
a2 


80 





ih eel RODUCTION 


Although numerical methods such as Runge-Kutta have been 
widely used to solve nonlinear differential equations by means 
of digital computers, the engineer does not always have access 
BPomomGaMpucet. Practical technigues to assist the engineer in 
obtaining approximate solutions are needed. 

Pipes (Ref. 1) developed a method of solving nonlinear 
erererential equations by a reversion technigue quite akin 
Bemaltgecbraic Series reversion. The technique uses Laplace 
transforms extensively in obtaining the approximate solution. 
Baycura and Brady (Refs. 2 and 3) have proposed and examined 
nonlinear Laplace transforms to obtain solutions to nonlinear 
differential equations. Tou, Doetsch, Pipes and Nowacki 
(Refs. 4, 5, 6 and 7) have proposed an iterative method using 
Laplace transforms in each iteration. Flake (Ref. 8) has pro- 
posed a cumbersome method uSing Volterra series and special 
transforms. 

The techniques were investigated and applied to driven 
second order systems with initial conditions. The systems 
were restricted to known stable systems, and with driving 
functions far from resonance. The method of Flake was found 
to be so impractical that it is only mentioned here for 


reference. 





ieee HOD Ss 


The various methods of obtaining approximate solutions to 


nonlinear differential equations were investigated. 


A. PIPES' METHOD 
Pipes has postulated a general nonlinear differential 


equM@ation OL the form 


, : 
OX Oe Foleo. Cena ae 0.” = k ew) —— 


where t is the independent variable, x is the dependent vari- 
able, k is a constant, 9(t) is a driving function, and the a. 
See runceLons Of the differential operator D, where 

eed cite (227k 2) 


Consider the special case of equation (2.1.1) where 





-> b; « T= Ot een (Zee) 
+) } e 

d= ia 

and the a. Aipemeenstanc Coerticients for i - 2, 3, «es, Nn. 


Equation (2.1.1) becomes 


(bLo+b,d ae thd REE 1) x ape ae 
dt” 
sdnx =k £0 


(2.1.4) 


Sreuscing the Operator notation 


Cea D: -seroeD \< + aX +... 


(2.1.5) 
tank zk ge 


Following Pipes, assume a series solution of the form 


x(t) = Ak +Awk +Aryk' see. (2.1.6) 





By substituting equation (2.1.6) into (2.1.1) and equating 
coefficients of equal powers of k, the various A, (t) terms 


may be found. They are 


A ® = gu) /a, (ene) 
Awe - a, A W/o, (DB 
AW een) za AWA, * aXe | (2.1.9) 


A,® =- Na. (oi 2hoA)s sa,NOAMsa, Nw | (Qaleet0) 


A.w-- \\za,\ MAO ASAD) 3a, A@A,03 NOh,®) 
+4a,N OAM s0,.N@ | 


(255 eels) 


Pipes (Ref. 1) lists additional terms. However, since a 
second or third approximation to the solution will generally 
be all that is desired, only the first five terms are given 
here, 

Be tairse order approximation can be found by solving the 
linear differential equation (2.1.7). Since (2.1.7) is a 
linear differential equation, it can be solved by the usual 
Eecsimiques Or Integration Or Laplace transforms, applying 
initial conditions. The remaining A, (t) are solved in like 


Manne tyexeeor the anitial conditions are zero. 


B. BRADY-BAYCURA METHOD 
Brady and Baycura have shown that for exponential 
functions and for small values of time the following approxi- 


mate expression holds 


L Ix'wr| x Ss XS) (21) 





Since the Laplace transform of 6(t) is known, a11 trans- 
forms form the terms of equation (2.1.3) are known ai.d thus 


(tes) can be solved for X(s). 


fee POU-DOEBTSCH-PIPES—-NOWACKI METHOD 

This method probably originated with Pipes, but was 
arrived at independently by the others. 

Assume the nonlinear system differential equation of the 


form 


iD a, D4 a +a. See ra. {ys Iya gee) a 4% (2205 el) 


Taking the Laplace transforms of both sides of equation 


(2.3.1) yields 
(9 a a.s"- Qu eG Vy Ke, yf = GG) C2322) 


where Q(s) is the Set of initial conditions. 


Mime nicse appLoximation assume the nonlinear term 


Ley. 4 a [=O , then 


Ys) = oS . (De os) 
) 
Siow: 
=o 
The inverse Laplace transform of Y¥, (s) yields the time 


Gena riesolutron Of the first approximation. 


The nonlinear term is included in the second approximation. 

Y 7 Qea1 Gor LH y hs 

Ae 
De ace 


kzo 


(25334) 





where f(t, Bae meperc a £Uncti0on of the first approxi- 
mation of y, (t). 
Successive approximations are found by continuing the 


same procedure. 





Toe APE ICA TIONS 


Nonlinear systems investigated were limited to second 
order with initial conditions and extended to systems with 


mercing functions. 


A. NONLINEAR SPRING 


Soveta-eea cimple Spring-mass system described by 


———— 


m d@x , Fay=O | (oo) 
4+ 


a 


The system is undamped, with an initial displacement of 
tO) = Xo and initial velocity set to zero. F(x) is the 


weenie torce, and for Ehis example, let its characteristic 


be 
- Ge = Lx rare (322) 
Letting all constants be unity, the system equation is 
nonlinear 
d* | 
ex x = OO (35 less 
a {> 
KO =X, 
: (eo eles 
xX) = O 


als Solution by Pipes Reversion Method 
Using the method outlined in Chapter II, Section A, 


equation (3.1.4) can be rewritten as 


10 





pet). 4 x 


= © 
(321.5) 
The as coefficients are seen to be 
2 

en =o Dy 

Cis, = | (32-36) 

Rk = | 

BY) =O 


The individual A. coefficients are then computed using the 


formulas (2.1.8), 


tet.) and (2,1,10): 


A @&= d@)/o, 


CAaler) 
2 
=O) 
Oa)A, Cole) 
Taking the Laplace transform of equation (3.1.3) and 


Mmeeessing initial conditions yields: 


A ®= <7 + 
! Saal (3.1.8) 
This well known transform gives 
A G@y= x,cost (Bee) 
As expected, the system is oscillatory. 
The next term A. (t) isum@ound frem equation (2.1.8) 
iN <3) ~&a A’ a Xe cost 
é Q, TY +) (GC yea ES) 


Aj (t) and all subsequent A, (t) are determined by the con- 


dition that they have zero initial conditions. 


tel 





Us sence IGemneity 
cos’ 4 = LCi sceszt) (3.1.11) 


and taking Laplace transforms of both sides of equation 





ie. 1.10) 
ee eee 
A = 2 {alsa Gora Coe * A) | (3.1.12) 
Born terms can be expanded into partial fractions: 
ole v \ 
page. - = - Pees es Yes (3.1.13) 
2 a = st ett si +4 ; 
tietng inverse Laplace transforms: 
x? 
A_@ = =e \\ - Zesst -lLcos2zt (3.1.14) 
7 2 S 3 


A, (t) is found in a Similar manner, and the solution is 


given by 
Xi) = Aw + Aw 4, Aw bes ae (2.1.6) 
Xe) e& XCost XL 1-2 cosk-b casat | 





Mieco (seit >) at is observed that the first 


secular? term occurs in the third approximation. 





25k 
secular terms are those involving t sin t, t cos fee Sil cy 


pit e LC. 


a2 





2. Solution by Brady-Baycura Transform Method 
tieriomEite wl lace Eranstorm Of equation (3.1.3) and 


applying initial conditions gives 
: 
S Xs) -8X, + XQ) + SXG) =O (3.1.16) 


In the usual manner, egquation (3.1.16) is Solved as a 


quadratic. After rearranging terms 





ae i ale AX,S ye | 
XG) oa | yy (1 + Gar (275 ib-1.7,) 


Neglecting the negative radical as physically impossible, the 


positive radical can be expanded as a binomial 





x 328 
eee eX, CS 


FH) SCRE NS 
Co Songs ice Pens) 


_ 


~—\ 
The first term is just x, cos t; however a Len | 


is not one of the common transforms. Reference 10 lists the 


following transform: 


L| ; i 2 eae a Y J n- 34 @.4) 


(ee _) 2 
(s*+ 0) ha /. ens one a (23 2re eo) 


where JIn-3/2 (wot) is the half order Bessel function. 
Further it is known from the theory of Laplace trans— 


forms that 


13 





-} n-!' 
ga on Cony = y= Ye 
[sq = w+ 7, eS We 2a a. 3 
=. ) | (3.1.20) 


Since £*) (9) Pe mean Oy). 2 2. ., n-l, the second term 
can be found: 


letting “=| 


) 
x, 4) = GAe)*( os = cost ) 





(32.203) 


Paemateem Caurying OUL the differentiations of equation 


(3.1.20), the second term of equation (3.1.18) becomes 





ass 
=x \atsint + teost | (3.1.22) 


SS 


Repeated application of the method will yield the 


marca term. 


The approximate solution to three terms is the truncated 


series 


KW) ® xX cost — Xe | Se sint + Peast | 
& 


ye\s wn ae Ss 3 ny 4 | an 
1% | S(ts t Peost) +S t Sys} +L cost 4 (3.1.23) 


14 





solute ton Dy fou-Doetsch-Pipes-Nowacki Method 


Putting equation (3.1.3) in the form of 025537) 


of Z 
CR eae xe = OC 


(eta 
where 
oe 
QA, 79 
GQ, =o 
on 
VQy 4.35 x (3.1.25) 
4 =o 
Eaiingeuaplace transforms of both sides 
=e 
Ss " _ 
MS, as -Qar Lh G44). > GG) (Sila 26) 
rb 
2 i L = 
Fe) eri) sxal |x | QO (3.1.27) 
Selving for X, (s) ase equation (2.322) 
~—OAesGO xs 
AS : — he Z 4 (S25) 
2, OAs 
hed 
Then taking inverse Laplace transforms yields the first 
approximate solution: 
xX We= % cost (eis529)) 


> 





The second approximation is found By applying equation 


fe 3.3) 
5} Q® + Gw ~ LEG, 0} 


aa 
2. as (Qala) 





AY = = ~ ay Poco} | (eee) 


Using the identity of equation (3.1.11), equation 


(3.1.31) becomes 


Si 
Xi =" 4 | -e5\4 re (3.1.32) 


Taking the inverse Laplace transform gives the second approxi- 


mation as 


X%, © =X Cost - XE i ae J coszt | 
As ie (312633) 


Repeating the method gives x, (t) and the result is the same 
acmequati1on (3.1.15). 

While this method appears to be different than Pipes 
method of Chapter II, Section A, consideration of the manipu- 
lation of the solution shows the two methods to be equivalent. 


Since the method of Chapter II, Section A is less cumbersom, 


it will be used hereafter. 


16 








4. Comparing the Methods 

Figures 1 and 2 are phase portraits of the solution 
as given by the IBM 360 computer. Runge-Kutta integration 
technique was used to obtain the solution. From these plots, 
it is clear that the system is oscillatory for x & 0.5 units 
and unstable otherwise. 

Figures 3, 4 and 5 compare the solutions by Pipes 
method and the Brady-Baycura technique for various initial 
Sonal tions. 

Themeectsno aemonstrate the relative accuracy of the 
approximations. The Brady-Baycura technic gives a better 
approximation in each instance; however, Pipes method is as 
Meecweatce fon Small initial conditions. The difference is 
Sawseq) by tne) Location of the secular terms in the solution. 
For both methods, however, the approximations are of little 
value for times greater than one unit. 

No attempt has been made to suppress the secular terms 
Which cause the approximate solution to diverge. It 1S gen- 
erally known that subharmonic resonance terms (submultiples 
of the natural frequency) exist in these nonlinear systems; 
they do not occur in this solution. The higher frequency Jump 


resonance terms do occur. 


B. NONLINEAR SPRING WITH DAMPENING 

The system of example 3.1 was modified to inelude Letrec1 ona. 
dampening. The resulting system equation takes on a velocity 
term. Letting all coefficients be unity, except the velocity 


term, the system equation becomes 


17 
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Zale 
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oe LOO 
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0.000 
OES 0 
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x = 
O 


RUNGE Kutta 


1.00000 
0.99750 
0.99002 
0,97 7E? 
0.96039 
0.76476 
0.20173 
-0.43604 


0.50000 
0.49906 


0.49626 


0.46675 
0.41002 
O747967 
=O sos 2 


OF 251000 
0.24961 
0.24844 
Oo21213 
0.1 lGes 


x = 
O 


x = 
O 
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TABLE I 


1 
Brady Pipes 
1.00000 LOCOCO 
0.99750 0.99854 
0.99002 0.99418 
O497762 0.98697 
0.96038 0.96424 
0.76469 0.86483 
0.22176 0.56088 
-0.23227 0.24276 
O50 
0.50000 0.50000 
0.49906 OS BED SNES 
0.49626 0.49678 
0.46674 Ona ito 
0.41002 0.42252 
0.18245 Oe22407 
=O, 075097 -0.02327 
O25 
0.25000 0.25000 
0.24961 0.24963 
0.24844 0.24850 
Ome 213 0.21369 
0.11214 Oeil? ES 





ax aX u 
oe eee = OC) 
VC (3.2.1) 





x 
‘6 
7 

tt 
% 
9 


AO) = © 
(3,2..2) 


The investigation included solving the system equation by 


Pipes and Brady-Baycura technigues. The final solution is 


found by letting the coefficient C approach zero in both 


methods. 


ae Solution by Pipes Method 


Usimaq operator notation, the a's AS 


Q, = Soccer e| 
Ql = \ 
B= 0 , (3.2.3) 


By letting Cc = 24, the algebra is simplified and the same 


result will obtain by letting 24 go to zero. 


Proceeding with the solution A, Oey a, and upon 


substituting and rearranging 


(3254) 


CDeeDnjAO=D 
Taking Laplace transforms and impressing initial COnadVeELons, 
equation (3.2.4) becomes 


24 





iS eS + \) A ts) ees ey 


(3-255 } 
or 
AY = eee Se 
Gaps siege (3.2.6) 


hoeoeence lt, formula 119, gives the inverse trans- 


form of equation (3.2.6) and 


A@- eS AYer) exp Ca ek : 





= (320) 
—CB-Ve Jexp bA- feu} | 
letting 28 eS to zero 
A a= x, exe Gd +exp Go | (3.2.8) 
re 
AwW= x.cost 
_ (3.2760) 


~The solution for A atter va modicum of algebra 


po 


becomes 


AW = = |. ~ 2% east 2, cos24 | (3.239) 


and substituting A. (t) ito ecuatiton (251.5) , A, (t) becomes 


2 a 
me — 

Kk) = —2— Joost ~Lapnot -Leostens2t | 
AS D+cDin i _ (3.2.10) 
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or 


L a Cret - A-CesZ2t - J cosat | 


Ao Sam 


3 (32.2501) 


Taking Laplace transforms, equation (3.2.11) becomes 


“Ss 


(Q - 3 \- fees 10 7 pease ap ara 
AO- % SS ZBS4) (S(O 2654) ) 
Sea eZ 


ome Ya o 16 
(S+4VS476s+ ) (S+N(S , 28+ D 


All inverse transforms of the terms of equation (3.2.12), 


Peeecmoe tne Second, can be found in Ref. 10. The second term 


is evaluated by convolution as follows: 


Let 
= 
pane 
S55 ( 3392 7 153)) 
CoE ens | 
2 S 428541 (See) 


then Zz, (t) meGOoS tt and 2. (t) = sin t when 2 goes to zero. 


Convolving 24 and 25 
E 
omc 2.) -| algae, (t=) dt. 


ae 
am a e 2 5 


The approximate Solution to three terms by Pipes 


becwmiguie 1S 
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eS ee ; 
XC) x xX,cost -— xX, \ == enst-\ cos Se | 
au % 3 


' | 
*S[S4 Sil -3 ¥ 2S emst ~ L Coszt-L cost} 4 Be 3a2nG) 
Sh 48 3 \ be 


tits Observed that the first secular term occurs in 
the third approximation. 
Beeeecolucion by Brady-Baycura Method 
The solution proceeds by taking the Laplace trans- 
form of equation (3.2.1), and impressing the initial con- 
ditions. The nonlinear Laplace transform is used for the 


x(t) term. 


SX (Ss) ¢ (Sy St MQ =X. (54 ¢) (3525 27 ) 


ByecevcGsLon Onc (Sexures technique, the left side of 
equation (3.2.17) can be made into a series by using the 
formulas of Ref. 9. 


Let equation (3.2.17) have the form 


7 CON bxXot bX 4... = ‘| (3m 2013) 


where 


be S4+eya 
a C32 tO) 


if 
vf 


Za 





then by reversion of the series 
pas) = Ty +t + Bayt ee 


where 


4 = XK Cs+e) 





e L3 Se 


ise pr GES) Cees ae 


Seeletting Cc = 2 
x (sS= esr ¢) 2 Se(sto) 
Sr 2AQs \ Se Z/s*1)° 


ws (s rey” Z 
(s+ 28s 1S 


A ‘ 


+ 


(3725.20) 
(saz) 
(3.3222) 
(3.72223) 
(3o224) 
(3622.25) 


The first term of equation (3.2.25) is the same as 


equation (3.2.6), and is just Xo Cos) t, 


SoluceiLOn. 


as in the Pipes 


Evaluation of the second term must be found by con- 


volution: 


Let 


28 





(322. 20) 


ZC) : <n +26% \ 
C3e25 27 } 


ZAI = = 245 
ay eral 


teaming inverse transforms, and letting “8 go to zero, 
(322. 23) 


Sa ney = Soa 


Car 229) 


Cy ya SOs Lt 


Then the inverse of the second term is the con- 


volution 
(3212 ..30) 


Zw =Z@) #20 ¥ B®) 
Again after a modicum of calculus, algebra and trigonometric 


substitutions, equation (3.2.30) reduces to just 
(Siz) 


Ale) : tS | = Cost rtemt | 


and to two terms, the Brady-Baycura technique yields 
(Se2eac) 


X(D ~ X,cost - = cook rtsint |. a 
eae 
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The technigue has lost its desirability even at this 
Peieesince tne Inverse transforms are so difficult and time 
consuming to evaluate. 

Seecompaning the Solutions 

The phase portraits given in me Sea, 7a and 7b 
show the apparent warping of the phase plane due to the 
inclusion of the velocity term. The instability of the 
system is seen to depend on both the value of C and Ky 

The Runge-Kutta solution is compared with the approxi- 
mate solutions of equations (3.2.16) and (3.2.32) in Figures 
8, 9 and 10. Table II lists selected values of the Runge- 
Kutta solution and the approximations. 

he iS Clear from the table that neither approximations 
are good even at 0.5 time units. It is further obvious that 
for values of the coefficient CL41.0, the system equation 
(3.2.1) will approach equation (3.1.3). Under the condition 
Saweeciatl,  cie techniques might have some value in approxi- 


Mating the transient state for small times (milliseconds). 


C. NONLINEAR SPRING WITH CUBIC TERM 


Consider again equation (3.1.1), but let 


Fas kes be (ue) 


Also let all coefficients be unity and the system have the 
same initial conditions as before. The nonlinear expression 


describing the system is 


ad x .% 4 K=O 


4 (33 220) 
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Time 
0.000 
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QO. 200 
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0.150 


weZ.00 


0.000 
0.050 
QO.100 
Oe 150 


0.200 


TABLE II 


x 


O 


Runge-Kutta 


#00000 
0.99754 
OF 991035 
Oo 700 0 


0.96290 


0.50000 


0.49908 


0.49638 


0.49200 


0.48604 


0.25000 
0.24962 
0.28849 
0.24666 


0.24417 


a 


x 


O 


O 


enw) 


Brady 
1.00000 


299719 
02988 77 
0.97478 


O35 95550 


ma O 0 


0.50000 
0.49898 
0.49594 
0.49089 


0.48384 


= O25 


Se 


0.25000 
OaZ4959 
0.24836 
0.24632 


0.24347 


Pipes 
1.00000 
0.99667 
0.98669 
0.97014 


0.94711 


0.50000 
0.49896 
0.49584 
0.49066 


0.48344 


0.25000 
0.24960 
0.24839 
0.24638 


0.24357 





ee 


Kis) = O (353) 


1. Solution by Pipes Technique 
Using operator notation, the as coefficients are 


identified as 
Q,= 41 
Q.= 9d | 
ae : (3.3.4) 
ly = \ 


B= oO 


The A, (t) coefficients are found next. From equation (2.1.7) 


AO= BoA, (reer) 


or 


(3752/5) 


i. ec 


Taking Laplace transforms and impressing the initial 
conditions on A, (t) the solution is immediately seen to be 


the same as that of the quadratic case 


A &= X%,cost (3.3.6) 


Substitution of as into equation (2.1.4), gives coefficient 


A, as zero. Ay is found from equation (2.1.4): 
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A ie) A\s. AQKQ) a) aS | 


(2ealare ) 





\ i, a Re 
= |> = «| (Sese0) 


A.& 3 


and after substituting for eee t, equation (3.3.7) becomes 


% =_— . 
Ke < 
Aa =| seek cee | (3.3.8) 


a “Dt+ | 


Taking Laplace transforms of both Sides 


A is) = ae eee —— (3.3.9) 
ae Se 


The inverse transforms for the factars of equation (3.3.9) 


are found in reference 11, and the solution of A, (t) is 


“> | 
AW: he | at snk s Leoel- Lecest | (3S LO) 


Since all even a. and thus even A. are zero, A, is 


zero. The coefficient A, after substitution, reduces to 


A\ - G\= - + | 20, MwA iw | 


(35, 95-135) 


3 





Pose palepiyeroreeie COCTTICients and manipulating trigono- 


metric identities, 


eit) = eens * Btsmt + Leost 
32041) 


— ese = Lcosst | 
4- 3 


and taking Laplace transforms 


A= Bx | er a, YS 


cen) is4Qyr (ey \) (s™ \ 5 (sta Lye 





-—- Yao Vet S 
(SF +9) es) (27425) (S405 


(cece) 


(Sos) 


The last three terms have readily available inverse trans- 


forms; however, the first and second terms require special 


Becnmniques. 
First term is evaluated by convolution: 


Let 


: =t} = 
ZA. lt) = - | =~ |- Cos gt 





Bt svat -s\o3t 
a 





1 \ 
aw- & ee. |* 


then 
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(225. sa) 


(325.135) 





1B| Ze) KZ) |: % [exax Bost) 3c] de 


(323 13¢)} 


= 4 Cost-Z25 cost -_34 sink 


tieemevaltatiton Of Ehe Second term is found in Appendix A: 
S ee 
Lo late ee: | Lernr Aaa 4 (3.3.13e) 


The remaining terms of (3.3.13) are found in reference 10. 
The result for A, (t) becomes 
= 
A ts) = 3“ i, 23 eost - 24¢c0s3t + onss +} 
Alo 


Atolgete Sins ls Se | 
a ) pase (nse) 


The approximate solution by Pipes method is the truncated 


series 


i 
Xe) we XCost - Xe Jat sint-+s cost .eos3t | 





<i ae oo. racosst s aosstl mesh 


+ A(sink Bom ae 34E east |: 
za 


4) 





Peeeomsect £rom equation (3.3.15) that the first 
secular term occurs in the second approximation. 
Peotone oy Beady-Baycura Method 
Upon taking Laplace transforms and uSing the non- 
linear expression of equation (2.2.2) and impressing the 


initial conditions, equation (3.3.2) becomes 


SX 4 C4 SKS) 'a-- <- (eS 2e3ie Gs) 


Considering equation (3.3.16) to be a series, the method of 
BevetstOn Outlined in Chapter III, Section B, paragraph 2 


can be used here also to obtain ajl the terms: 





’ sé 4 . (35357 ay) 


i, 2 (3.3.17b) 


BW 
uv 
i 


(Sex a 


\ = 2 
‘bs (2b ee ane = oe (3273..17-6) 


The solution in terms of X(s) is 


n aye > 
ee td 
eat fen) (stevy? (roe) 


The inverse transform of the first term is just x, COS ts 
The second term must be solved using the technique outlined 
in Brady (Ref. 2) and in the example of Chapter III, Section 
A, paragraph 2. The method is straightforward, giving the 


second term as 
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— X%& Se} stewt - eat atest | 


a eee) 


The second term (3.3.19) yields the higher order secular 
terms te cos t and t? Sin t which will be small contributions 


mor t / 1. 


It is also known from equation (3.1.19) that 


g 


a | wit Jy 


Susy 2% Tuy. | (se3ac0) 


By Observing the result in second term, one sees that 
Swomeiigier Order (t®) secular terms will result from the 
third term of equation (3.3.18), and further that the 
coefficient for this third term is small (less than 


oe x oO); unerefore, little is gained by solving equation 


(en. 20) . 


The approximate solution by the Brady-Baycura 


Peecmnlague 1S 


PA eee OSG = es | 4ekewt - Pormt 4 Feast vv. (3.3.21) 


Eeeecotipakbicon Of Solutions 
Figures 11 and 12 show the phase plane solution for 
equation (3.3.2). Here the solution is seen to be stable for 


eile Xo° The system always remains in a limit cycie. 


43 





Pio cOnooserveq from Figures 13, 14 and 15 that 
Pomelange@ initial Conditions, the Brady-Baycura method gives 
Moonen dppLoximatiton £O the solution for early time inter- 
mes. For the case of X Shatner ioe Sometmiod Gives better 
results. This is caused by the location and number of 
Becttlar termS appearing in the solution, Which causes it 
memaLlverge. 

Selleeecd values tor the graphs of Figures 13, 14 and 


15 are given in Table III. 


D. NONLINEAR SPRING WITH DRIVING FUNCTION 

Consider the system of Chapter III, Section A, but now 
as a forced system. The system nonlinear differential 
equation is 


max , Fa= Beosdt 
Qt (3.4.1) 





where 


Fud- lex sb (3.4.2) 


Letting all coefficients be unity, equation (3.4.1) becomes 


au 
es 4% xg = Warns dt (Saree) 
KX bo) = Key 
X@)= Oo 
(3.4.4) 
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Runge-Kutta 
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Oro sO 17 
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On 23570 


0.50000 


0.49922 
0.49688 
0.48064 
0.42461 


Vie Zoo 7Z 


0.25000 
0.24967 
Oe23 0415 
On.b7 923 


0.12946 


x 


ORS I 207 


0 


OR Ae 


Bo) 


PABLE LID 


Pipes 
1.00000 
0.99768 
0.99073 
0.97920 
0.94277 
0278095 


OZ oo 


0.50000 
0.49922 
0.49690 
0.48077 
0.42498 


One 2707 


O25 000 
QO. 24967 
Oe Zs516 
O27 930 


0.12947 


OZ 009502 


Brady 
1.00000 


O-297 50 
O. 22002 
OF Tet 
02 2353) 
On o2s 


On 19357 


0250000 
0.49922 
0.49688 
0.48063 
0.42444 


O22 260% 


OZ O00 
0.24967 
0.235605 
O27. 
Dalz I6e 
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Pes olmeLom OY Pipes Method 
The a. Gecthicients are these of equation (3.1.6), 
except here it is seen that @(t) aS cos Wt. 


The solution proceeds as before 


= (+) = BeosOt 
A te) p ja. een (304 5)) 


Taking Laplace transforms and impressing initial conditions, 


AIS) = a aSee | (3.4.6) 


+ " wy ero 
s 


The terms of A, (s) have transforms in Ref. 10. 


A la)= cost + 2 cost casas | (3.4.7) 
ae. 
= ky oct - ke, Cost (3.4.8) 
where 
eee tS 
Se 
(3.4.8a) 
_ SB 
k, -  y 


The solutions to A. (t) and A, (t) are straighforward 
using the method outlined in Chapter II, Section A. Thus 


A K)= les PCs Ry eset + ko Cos (Qt 


+ kyo Sssn)t-k, eoszast | (3.4.9) 


DL 




















cena Aa | ae \ = en 
f a3. Lote > —2 4(4ut-.) 

ee ER (3.4.9a) 
O.(,-2) 

iene oe es 
DS, G,4+2) 

kh, ~ 
ACTED 


The solution for A, (t) is given in Appendix B. 


The complete solution by Pipes method is given by 


XW= Aw skis Aya ww, (3.4.10) 
2. Setuetoneoy —Blady-Baycura Technique 


Equation (3.4.3) can be transformed as in Chapter 
ml, 


Section A, paragraph 2 uSing the nonlinear transform 
iE@r eee 


After impressing the initial conditions, 


the 
transformed system equation becomes 
Q : 
Ee Coys a) Ko) = X Se 1s | (een) 
Sy Le 


Letting equation (3.4.11) be a series, then by the method of 


algebraic series reversion 


XO: A Cas + 2s VA. (x4 a8s Na ice 2)) 
Soe 


st 4 St 
d 
where the A. coefficients are 


a2 





a? eae 
As -_s 
2 (ait (Sa4013) 
A\. = Bice 
(says 
Thus 
Zz 
iMieaj= | ~ Bs -|--s eS Se 
exh SPO) (s°+1)3 S44 i 
(3.4.14) 


The first term of the series in equation (3.4.14) 


is easily evaluated by performing the algebra indicated and 
Let the first term 


Pormigecne transform pairs of Ref. 10. 
be FL (s), then 
= 
Fid- ee en ee eee : 
: eet (S414 (3.4.15) 


and 


= k cost - kosot ora a0 7) 


where 
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ae | (see) 


The result is not unlike the first approximation 
Socainead im the previous method by Pipes. 


Let the second term be Fi(s), then 


Bs ¥ jens = : | 
i Ss fo sree 





: | KS rn 2BxX,S° 
(s+ 1)* (Sr 1)* 4 QP) j a ice LTRD (354220) 


The terms of F.(s) are all uncommon transforms. The 
ZeecwrLgue OL finding the inverse transforms is as follows: 


Since 


° lea + tamt- Voost ) 





(3.4.21) 


As shown in Appendix A, the derivative method of finding 


inverse transforms is extended to equation (3.4.21). 


= fa tse -tene)| a Le\atant+Penct | (3.1. 21a) 


endethe inverse transform of the first term of equation 


(364520) is 








The second term of F,(s) is evaluated by convolution. Except 


mora Constant, it is the first term of £.(t) convolved with 
=y | \ 

keer 

or 


ZBx| YelBtent fenet |x| enagt | (3.4.23) 


The expression (3.4.23) can be evaluated as 


=. BS lea 4, CESS cCosmt +k wtsrmit 4 RE 








\G WB, 
~GL_ sn @X vt - eT cae VB, ye | (sade 
Ld = \ ATID 
where 
eee = St 
ee Sr) es 
a wrese , 4 Peat: 
aes ee Ae 
Gays7 So EP OG 
ee eee ee sy (3.4.24a) 
LO. | U>-' CLds+ ye [sy 
LJ Z Ret, oe 
Ee LX* | 
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Exeepietor a Constant, the last term of F.(s) is 
Simply the term of (3.4.24) convolved again with 1 Swat 
Pie cesule 1S Lound in Appendix C. . 

Clearly the method is too cumbersome to carry beyond 
this second approximation. The approximate solution by the 


Brady-Baycura technique is 


eyes 
% Kt) 7,0 48 was 4 (Sx4e26) 


Sa compaxring the Solutions 


The time plots of Figures 16, 17, 18 and 19 exhibit 
eae relative accuracy of the approximations to the Runge- 
Kutta solution, for various combinations of Kos B and GU). 

From Table IV it is apparent that the approximations 
meen good Only for the first initial time units. Pipes method 
gives only a slightly better approximation than the Brady- 
Pemseuna technique. For the case u) | - 5.0, it is observed 
Bicdtmele ApoLOximations by both methods improve. This is 
in keeping with the original restriction of the system to 


PeLeGiniag functions far from resonance. 


E. NONLINEAR SPRING WITH CUBIC TERM AND FORCING FUNCTION 
Assume the describing equation to be 
% 3 
cx er = BCosW 
ae (ease) 


Caeheinatcial conditions 
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Kw\= X, 
alte = O (54:5%,02) 


ale 


SOluEton Dy Pipes Method 


The solution proceeds in the usual way be identifying 


a.'s 
1 \ 
Q= Da\ 
Qo: o 
Q.~ \ 
e oe 


(3252) 
Die) = Beosws | 


Using equation (2.1.7), the first approximation Aj. is found 


Pyestaplace transform techniques. 


Atay = SS (3.57.5) 
\ ‘Sas \ 
8 (3.5.4) 
\ sh \ B41) Gk) 
The solution is the same as for the quadratic case 
A G)\= kcosk - , Coa dat (3.5.5) 
where 
Ca mee 
aN 
(3297/0) 
kB 
LS*— } 
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The coefficient A. (t) Was) SCen) prevroueiy to bewzero. 


Solving for A, (t) is straightforward and the result is 


Eat) = [tent a e Sr by , ae 
wi +) = 
_ Alqas ayy] 





—_ <a 
41 Quss-1] qaaya) as weary +o coset 


y Sey ky endork, BEL codto,r2¥ 
4\(o-7- | 4-[(to425- 1 | 


(325.7) 
BRR, cog Gut) meee cs Cos(2,-1)+ 
~ 41s Yo) Herd] 
a 
+ Se Re COS Sut 
Aloe) 1) 
where 
Ry? 3eY x % kha 
4 
(32526) 


Ras Hbrk + HR 


Miewsolution of the undriven system of Chapter III, Section 
C should demonstrate the Mincoacloliity Or Carrying ene 
approximation beyond A, (t). Furthermore, the first secular 


term in A, (t) will cause the solution to diverge. 


2. Demon oy Isrady-Baycura Technique 


Taking Laplace transforms of equation (3.5.1) and 
impressing the initial conditions, the system equation 


becomes 


oe 4 (S4 NXE) = ox, + —_ BS — (3.5.9) 


See 


63 





Using the algebraic reversion of series technique to solve 


for X(s), let the left side of equation (3.5.9) be a series 








ee se 
Ny or Or Wx© Oe ee (35510) 
then 
RGs Ey | ~ 3 
s) 4%) Sy By * Day's res (3. cael 
NX 
= ee 4 Ss > B (sor Bs \s Ss soar) 
aS seco 
Where from Ref. 9 
BR - RE. = 
FT », > +) (3. ae) 
ane 
> (Seay 
i. te Rb) -- Se (3.5.15) 
> b+ Gros 
Then to first order approximation 
X (s) = By (Bese ia) 
= =e SS ron 
and 
XG) = b cost ~k onsd-t (35218) 
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where 
: RB 


(3.5.18a) 
aS 


LQ \ 
AS expected, the result is equivalent to Pipes MeEnod for 
the first approximation. 


Since Bo is zero, the third term, F(s), of the 


series is found 


- 
= & ISs 


z (s“41yh Srey Gece, 
foes ee | Kee 
ee eee LC o ai 
n= 4. Se 
Gay ar) ($+.at) (Frat) (st 
. aos (35.20) 
(Sr ay)(ss yt 
pieeinverse transform of the first term as given in 
eamation (3.3.18). 
2h S \_ - 
vg __ = — ist sit «9 Lect -sut| 
Ge © | (G5e20) 


Except for the constants, the successive terms of F,(s) must 
be evaluated by successive convolution of = sind. t with 
Best igmt Side of equation (3.5.21). eae eneanain’ the 
evaluation of F.(s) will lead to terms of cos t, modulation 
terms and higher resonance terms, it will also contribute 
numerous secular terms, and thus cause the solution to di- 


verge more rapidly. Again the technique has lost its desir- 


Bomeitey oy providing cumbersome transforms. 
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Commemene £irSt term only of Fj(s), Ene Solution is 
auven by 
xU~ k, Cost -e,cosat - x3 piste +9 Penst 


4B ms 
~ trant|s.. ‘ (3.5.22) 


3. Comparing the Solutions 

feguees 20, 21, 22 and 23 are graphs of the approxi- 
mate solutions and the true solution for various combinations 
of 7 2 anda .. Table Vv lists selected YValueccmironeene 
Serucions. 

ioe Cpesvejas elsiurengkeis ethers that except aepreleite acme eta 
case of Figure 20 where Xo a epeandeay = = 1.5, that Pipes 
Metriod provides only a Slightly better approximation than 
the Brady-Baycura Method. This is due to the number of 
terms generated by the Pipes Method. Figure 21 shows that 
the condition of the system operating far from resonance, 
Pipes Method gives a better solution for a longer period of 


mine we bOth approximations improve as Xo and B become small. 


66 








Oo6O 








0.50 
PIPSS 

0.40 

©o.30 

ree) \ 

RUWGE=KULPA ——7" \ 
7 \ 
. E 050 % 

~O.40 


3 


FIGURE 20. Time solution of x + x * x” ~ 0.5 cos 1.5t, 
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FIGURE 21. Time solution of x t+ xt ce = 0.5 cos 5.0t, 


Xo = (039 
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FIGURE 22. Time solution to SG _ = 0.5 cos 1.5t, 
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FIGURE 23. Time solution to Seed 3c) >: ee = 0.25 cos 1.5t, 
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Runge-Kutta 
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Bye CONCLUSIONS . 


It has been demonstrated that these various methods of 
obtaining approximate solutions to nonlinear systems with 
bela Conditions are accurate only during the early time 
Pmieemvals. Further, it was shown that this degree of accuracy 
May be Satisfactory for engineering applications. 

It has also been shown that the Brady-Baycura technique 
of nonlinear transforms usually provides a more accurate 
solution; however, the ease and facility of Pipes Method 
might generally make this a more desirable method. 

The method of handling so-called secular terms is still 
unresolved. .It is suggested that they may be eliminated by 
fetting the frequency also be a series such that the con- 
stants of the series are chosen to eliminate resonance con- 
ditions and thus eliminate the secular terms. The method 
would be similar to that given by Pipes in Ref. 12. 

Both methods show that the higher order harmonic terms 
exist in the solution; however, the subharmonic terms fail 
to occur in the first few terms of the approximations. 

Since the problem is an initial value problem, it appears 
that the solution could be improved by stopping the problem 
solution before it diverges from the true solution (t .5), 
then restarting with a new set of initial conditions. The 
effects of x(o) would have to be investigated. The value of 
such a method seems dubious since Runge-Kutta analysis on the 


digital computer seems an easier method. 
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APPENDIX A 


Ss 
+ 1) 


INVERSE LAPLACE TRANSFORM OF 


(s? 


3 


The derivative technique is used: 


From Ref. LO, 


til . L\sut~ keost | 





also 





or 


ae -i4d\ ! ae 
(sty...) 4ds | Qty 


since 


P-L a(n) we PT 
. E ‘|= ea) ee eps 


where n is the aad derivative. Thencombining equations 


(A-4) and (A-3) 


Clete | 3 (dies 
Set te | 
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ae 


(A-2) 


(A-3) 


ae) 


(A=) 


(A-6) 





APPENDIX B 


EVALUATION OF Ay IN FORCED QUADRATIC SYSTEM 


From equation (2.1.5) A, becomes 


A et, CN " A i) ae) 


(Baas) 
oS 
Perer substituting for ay A, and A. and combining terms; 
A Be) ware Qest — ie. fi eons - le gCBS9,t 
WF te, goss rk a bens 2k + ky cost 
+ Rey ers “yt eS & COS (244k SONS oUt 
. (BSc) 
= b, eos Bot i ba keg | 
where 
lon = 2k, - Be a bee i b key 
ls 
4 = 2k,k,- Daze a bk, ~ bby 
pe 
kis = R ke ‘ RR (B-3) 


lo 


kn = Bee t RR, 
G 
en = k le, Bs a ke 


Riko +k ky 
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Except for the first term of Ax, the transforms are all of 




















the form: 
= 
(say) (Sto) (B-4) 
For which the inverse transform is known: 
ai |erst -cosat | 
Ls) . (B25) 
The first term has the transform: 
a 
B-6 
which has the inverse 
& sint (B-7) 
ee 
Using expressions (B-5) and (B-7), A, becomes 
~ 4 
AXz\t)= | Betant ~( fates R. i hg - eae ~ key 
S38 QD Goh &eh 
ne ki kg = ‘ev leg ~~ Bow + es =< kek, Cost 
Ge Gea) Ge Gah 
> 
+ Bkeg costt + RY nest + ka onsat 
3 4-2, Lt 
(B-8) 


—s ces(w-2)+ _ ky ons(otz)t- Reiley casos Jt 
Cor 2) = (a2j-1 Gx) 


a ler ky ens (wart + ea ans@ls tt a @ps.(2u -1Yt - bk, Mee sone 
(0-2) Aly(L9,+)) Aw (v1) Awa 
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APPENDIA C 


SOLUTION OF LAST TERM OF F.(s) BY BRADY-BAYCURA METHOD 
Eeteenierexr pression in parenthesis of (3.4.24) equal Z, (t) 


and + sin Peco at), then 


O 2 | 
Oo 


Zz (t) = fw rt G)4.., (C-1) 


thus 


2 





\2 Hz) - $ORZW +4,.W KZM¢... 


32 W; 





(Ca) 


aor weconvOluciOns make up the last term. Since each of 
these convolutions is of the type previously evaluated, the 


solution is now straightforward and the last term becomes 


& is Rg E> £4), cost , _\) tos@oyt _ 20, Cos O 


us* —| Se. \ ast 


* bas re ten - Cas SO,t ) — ae % CoSs3s0 + | 


- any i \ = Cost 4 \ es 
z ee 3 Baw) * arp Gg) oO, 


+k, ZW, Cost Peete oe 
a ae) us —| >t 








Sin (2 4r)t 














4 sin@e~)t + cos (Zak +_) cosG@ay DY 
to) Lo \ ud,-\ 
¥4tf_\ 1) Vat 4 2 [sind - sina} 
(ud Deas (ca, ¥1\% 
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= 





~ Cost — cos(20, +1)t Pe wA + Z.cosiae| 


es Ssiat x Sin (ZS Nt ~ Cost + Cos@wo-)t ~ZSin OL 
Coe 3h 








a: 
Beco as)| oe \( ae — ape 


+ ps. a Saar kL An aay. 4 esata WEI Cos Ok 
3 wo) L,-) ud,- Do (23 = Vs 


~ 3 eee ast ef SN \enyst (c-3) 
oy [ sds1 W.-1 Qs) Bas 
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